Taking Multiple Infections of Cells and Recombination into Account Leads to Small Within-Host Effective-Population-Size Estimates of HIV-1

Whether HIV-1 evolution in infected individuals is dominated by deterministic or stochastic effects remains unclear because current estimates of the effective population size of HIV-1 in vivo, N-e, are widely varying. Models assuming HIV-1 evolution to be neutral estimate N-e similar to 10(2)-10(4), smaller than the inverse mutation rate of HIV-1 (similar to 10(5)), implying the predominance of stochastic forces. In contrast, a model that includes selection estimates N-e>10(5), suggesting that deterministic forces would hold sway. The consequent uncertainty in the nature of HIV-1 evolution compromises our ability to describe disease progression and outcomes of therapy. We perform detailed bit-string simulations of viral evolution that consider large genome lengths and incorporate the key evolutionary processes underlying the genomic diversification of HIV-1 in infected individuals, namely, mutation, multiple infections of cells, recombination, selection, and epistatic interactions between multiple loci. Our simulations describe quantitatively the evolution of HIV-1 diversity and divergence in patients. From comparisons of our simulations with patient data, we estimate N-e similar to 10(3)-10(4), implying predominantly stochastic evolution. Interestingly, we find that N-e and the viral generation time are correlated with the disease progression time, presenting a route to a priori prediction of disease progression in patients. Further, we show that the previous estimate of N-e>10(5) reduces as the frequencies of multiple infections of cells and recombination assumed increase. Our simulations with N-e similar to 10(3)-10(4) may be employed to estimate markers of disease progression and outcomes of therapy that depend on the evolution of viral diversity and divergence.

Stochastic Simulations Suggest that HIV-1 Survives Close to Its Error Threshold

The use of mutagenic drugs to drive HIV-1 past its error threshold presents a novel intervention strategy, as suggested by the quasispecies theory, that may be less susceptible to failure via viral mutation-induced emergence of drug resistance than current strategies. The error threshold of HIV-1, mu(c), however, is not known. Application of the quasispecies theory to determine mu(c) poses significant challenges: Whereas the quasispecies theory considers the asexual reproduction of an infinitely large population of haploid individuals, HIV-1 is diploid, undergoes recombination, and is estimated to have a small effective population size in vivo. We performed population genetics-based stochastic simulations of the within-host evolution of HIV-1 and estimated the structure of the HIV-1 quasispecies and mu(c). We found that with small mutation rates, the quasispecies was dominated by genomes with few mutations. Upon increasing the mutation rate, a sharp error catastrophe occurred where the quasispecies became delocalized in sequence space. Using parameter values that quantitatively captured data of viral diversification in HIV-1 patients, we estimated mu(c) to be 7 x 10(-5) -1 x 10(-4) substitutions/site/replication, similar to 2-6 fold higher than the natural mutation rate of HIV-1, suggesting that HIV-1 survives close to its error threshold and may be readily susceptible to mutagenic drugs. The latter estimate was weakly dependent on the within-host effective population size of HIV-1. With large population sizes and in the absence of recombination, our simulations converged to the quasispecies theory, bridging the gap between quasispecies theory and population genetics-based approaches to describing HIV-1 evolution. Further, mu(c) increased with the recombination rate, rendering HIV-1 less susceptible to error catastrophe, thus elucidating an added benefit of recombination to HIV-1. Our estimate of mu(c) may serve as a quantitative guideline for the use of mutagenic drugs against HIV-1.

Estimating the fraction of progeny virions that must incorporate APOBEC3G for suppression of productive HIV-1 infection

The contest between the host factor APOBEC3G (A3G) and the HIV-1 protein Vif presents an attractive target of intervention. The extent to which the A3G-Vif interaction must be suppressed to tilt the balance in favor of A3G remains unknown. We employed stochastic simulations and mathematical modeling of the within-host dynamics and evolution of HIV-1 to estimate the fraction of progeny virions that must incorporate A3G to render productive infection unsustainable. Using three different approaches, we found consistently that a transition from sustained infection to suppression of productive infection occurred when the latter fraction exceeded similar to 0.8. The transition was triggered by A3G-induced hypermutations that led to premature stop codons compromising viral production and was consistent with driving the basic reproductive number, R-o, below unity. The fraction identified may serve as a quantitative guideline for strategies targeting the A3G-Vif axis. (C) 2013 Elsevier Inc. All rights reserved.

Scaling law characterizing the dynamics of the transition of HIV-1 to error catastrophe

Increasing the mutation rate, mu, of viruses above a threshold, mu(c), has been predicted to trigger a catastrophic loss of viral genetic information and is being explored as a novel intervention strategy. Here, we examine the dynamics of this transition using stochastic simulations mimicking within-host HIV-1 evolution. We find a scaling law governing the characteristic time of the transition: tau approximate to 0.6/(mu - mu(c)). The law is robust to variations in underlying evolutionary forces and presents guidelines for treatment of HIV-1 infection with mutagens. We estimate that many years of treatment would be required before HIV-1 can suffer an error catastrophe.

Rule-based modelling of iron homeostasis in tuberculosis

To establish itself within the host system, Mycobacterium tuberculosis (Mtb) has formulated various means of attacking the host system. One such crucial strategy is the exploitation of the iron resources of the host system. Obtaining and maintaining the required concentration of iron becomes a matter of contest between the host and the pathogen, both trying to achieve this through complex molecular networks. The extent of complexity makes it important to obtain a systems perspective of the interplay between the host and the pathogen with respect to iron homeostasis. We have reconstructed a systems model comprising 92 components and 85 protein-protein or protein-metabolite interactions, which have been captured as a set of 194 rules. Apart from the interactions, these rules also account for protein synthesis and decay, RBC circulation and bacterial production and death rates. We have used a rule-based modelling approach, Kappa, to simulate the system separately under infection and non-infection conditions. Various perturbations including knock-outs and dual perturbation were also carried out to monitor the behavioral change of important proteins and metabolites. From this, key components as well as the required controlling factors in the model that are critical for maintaining iron homeostasis were identified. The model is able to re-establish the importance of iron-dependent regulator (ideR) in Mtb and transferrin (Tf) in the host. Perturbations, where iron storage is increased, appear to enhance nutritional immunity and the analysis indicates how they can be harmful for the host. Instead, decreasing the rate of iron uptake by Tf may prove to be helpful. Simulation and perturbation studies help in identifying Tf as a possible drug target. Regulating the mycobactin (myB) concentration was also identified as a possible strategy to control bacterial growth. The simulations thus provide significant insight into iron homeostasis and also for identifying possible drug targets for tuberculosis.

How does contrasting dependence of impurity-atom diffusivity on the density of host disordered medium arise?

We report results of molecular dynamics investigations into neutral impurity diffusing within an amorphous solid as a function of the size of the diffusant and density of the host amorphous matrix. We find that self diffusivity exhibits an anomalous maximum as a function of the size of the impurity species. An analysis of properties of the impurity atom with maximum diffusivity shows that it is associated with lower mean square force, reduced backscattering of velocity autocorrelation function, near-exponential decay of the intermediate scattering function (as compared to stretched-exponential decay for other sizes of the impurity species) and lower activation energy. These results demonstrate the existence of size-dependent diffusivity maximum in disordered solids. Further, we show that the diffusivity maximum is observed at lower impurity diameters with increase in density. This is explained in terms of the Levitation parameter and the void structure of the amorphous solid. We demonstrate that these results imply contrasting dependence of self diffusivity (D) on the density of the amorphous matrix, p. D increases with p for small sizes of the impurity but shows an increase followed by a decrease for intermediate sizes of the impurity atom. For large sizes of the impurity atom, D decreases with increase in p. These contrasting dependence arises naturally from the existence of Levitation Effect.

Multiple molecular alterations in phosphodegrons 1-3 within AAV2 capsid demonstrates higher hepatic gene transfer efficiency

The success of AAV2 mediated hepatic gene transfer in human trials for diseases such as hemophilia has been hampered by a combination of low transduction efficiency and a robust immune response directed against these vectors. We have previously shown that AAV2 is targeted for destruction in the cytoplasm by the host-cellular kinase/ubiquitination/proteasomal degradation machinery and modification of the serine(S)/threonine(T) kinase and lysine(K) targets on AAV capsid is beneficial. Thus targeted single mutations of S/T>A(S489A, S498A, T251A) and K>R (K532R) improved the efficiency of gene transfer in vivo as compared to wild type (WT)-AAV2 vectors (∼6-14 fold). In the present study, we evaluated if combined alteration of the phosphodegrons (PD), which are the phosphorylation sites recognized as degradation signals by ubiquitin ligases, improves further the gene transfer efficiency. Thus, we generated four multiple mutant vectors (PD: 1+3, S489A+K532R, PD: 1+3, S489A+K532R together with T251 residue which did not lie in any of the phosphodegrons but had shown increased transduction efficiency compared to the WT-AAV2 vector (∼6 fold) and was also conserved in 9 out of 10 AAV serotypes (AAV 1 to 10), PD: 1+3, S489A+K532R+S498A and a fourth combination of PD: 3, K532R+T251. We then evaluated them in vitro and in vivo and compared their gene transfer efficiency with either the WT-AAV2 or the best single mutant S489A-AAV2 vector. The novel multiple mutations on the AAV2 capsid did not affect the overall vector packaging efficiency. All the multiple AAV2 mutants showed superior transduction efficiency in HeLa cells in vitro when compared to either the WT (62-72% Vs 21%) or the single mutant S489A (62-72% Vs 50%) AAV2 vectors as demonstrated by FACS analysis (Fig. 1A). On hepatic gene transfer with 5x10^10 vgs per animal in C57BL/6 mice, all the multiple mutants showed increased transgene expression compared to either the WT-AAV2 (∼15-23 fold) or the S489A single mutant vector (∼2-3 fold) (Fig.1B and C). These novel multiple mutant AAV2 vectors also showed higher vector copy number in murine hepatocytes 4 weeks post transduction, as compared to the WT-AAV2 (∼5-6 Vs 1.4 vector copies/diploid genome) and further higher when compared to the single mutant S489A(∼5-6 fold Vs 3.8 fold) (Fig.1D). Further ongoing studies will demonstrate the therapeutic benefit of one or more of the multiple mutants vectors in preclinical models of hemophilia.

Divvying up an incubator: How parasitic and mutualistic fig wasps use space within their nursery microcosm

Differential occupancy of space can lead to species coexistence. The fig-fig wasp pollination system hosts species-specific pollinating and parasitic wasps that develop within galls in a nursery comprising a closed inflorescence, the syconium. This microcosm affords excellent opportunities for investigating spatial partitioning since it harbours a closed community in which all wasp species are dependent on securing safe sites inside the syconium for their developing offspring while differing in life history, egg deposition strategies and oviposition times relative to nursery development. We determined ontogenetic changes in oviposition sites available to the seven-member fig wasp community of Ficus racemosa comprising pollinators, gallers and parasitoids. We used species distribution models (SDMs) for the first time at a microcosm scale to predict patterns of spatial occurrence of nursery occupants. SDMs gave high true-positive and low false-positive site occupancy rates for most occupants indicating species specificity in oviposition sites. The nursery microcosm itself changed with syconium development and sequential egg-laying by different wasp species. The number of sites occupied by offspring of the different wasp species was negatively related to the risk of syconium abortion by the plant host following oviposition. Since unpollinated syconia are usually aborted, parasitic wasps ovipositing into nurseries at the same time as the pollinator targeted many sites, suggesting response to lower risk of syconium abortion owing to reduced risk of pollination failure compared to those species ovipositing before pollination. Wasp life history and oviposition time relative to nursery development contributed to the co-existence of nursery occupants.

Parasites exert conflicting selection pressures to affect reproductive asynchrony of their host plant in an obligate pollination mutualism

1. Plant reproductive phenology is generally viewed as an individual's strategy to maximize gamete exchange and propagule dispersal and is often considered largely dependent on patterns of floral initiation. Reproductive phenology, however, can be affected by proximate responses to pollinators, parasites and herbivores which could influence floral longevity or fruit development time. 2. We examined the influence of insect interactants on within-plant reproductive phenology in the fig-fig wasp nursery pollination mutualism in Ficus racemosa (Moraceae). Most figs support a wasp community comprised of a mutualistic pollinator, with several host-plant-specific non-pollinating herbivorous gallers and parasitoids. These wasps reproduce within enclosed inflorescences called syconia, which develop into fruit after pollination. While different wasp species oviposit into syconia at varying times during its ontogeny, all wasp progeny are constrained to exit syconia simultaneously just prior to fruit ripening. Developing larvae of early-ovipositing wasps may hasten syconium ontogeny through formation of earlier and larger nutrient sinks, whereas larvae of late-arriving parasites may lengthen syconium ontogeny to complete their development successfully. Seeds are also important nutrient sinks. The number of seeds and the type and number of developing wasps may therefore be expected to influence syconium development times, thereby affecting the reproductive synchrony of syconia on a plant. 3. Observations on naturally pollinated and parasitized syconia indicated that their seed and wasp content affected syconium development time. Experimental manipulations of syconia to produce only seeds or various combinations of wasps confirmed this finding. Early-ovipositing galler progeny reduced syconium development times, while gallers ovipositing concurrently with pollinators had no effect on syconium development. Late-ovipositing parasitoid progeny, the presence of only seeds within the syconium, or delayed pollination increased syconium development time. The differential development of syconia, which was influenced by mutualistic or parasitic progeny, accordingly contributed to within-tree reproductive asynchrony. 4. Synthesis. Individual reproductive units in fig trees called syconia, which also function as brood sites for pollinating and parasitic fig wasps, have plastic development durations dependent on pollination timing and species of wasps developing within them. Syconium development times are a likely compromise between conflicting demands from developing seeds and different wasp species.

Beyond classical dynamic structural plasticity using mesh-free modelling techniques

The problem of modelling the transient response of an elastic-perfectly-plastic cantilever beam, carrying an impulsively loaded tip mass, is,often referred to as the Parkes cantilever problem 25]; The permanent deformation of a cantilever struck transversely at its tip, Proc. R. Soc. A., 288, pp. 462). This paradigm for classical modelling of projectile impact on structures is re-visited and updated using the mesh-free method, smoothed particle hydrodynamics (SPH). The purpose of this study is to investigate further the behaviour of cantilever beams subjected to projectile impact at its tip, by considering especially physically real effects such as plastic shearing close to the projectile, shear deformation, and the variation of the shear strain along the length and across the thickness of the beam. Finally, going beyond macroscopic structural plasticity, a strategy to incorporate physical discontinuity (due to crack formation) in SPH discretization is discussed and explored in the context of tip-severance of the cantilever beam. Consequently, the proposed scheme illustrates the potency for a more refined treatment of penetration mechanics, paramount in the exploration of structural response under ballistic loading. The objective is to contribute to formulating a computational modelling framework within which transient dynamic plasticity and even penetration/failure phenomena for a range of materials, structures and impact conditions can be explored ab initio, this being essential for arriving at suitable tools for the design of armour systems. (C) 2014 Elsevier Ltd. All rights reserved.

Identifying regions for conservation of sloth bears through occupancy modelling in north-eastern Karnataka, India

In the absence of information on species in decline with contracting ranges, management should emphasize remaining populations and protection of their habitats. Threatened by anthropogenic pressure including habitat degradation and loss, sloth bears (Melursus ursinus) in India have become limited in range, habitat, and population size. We identified ecological and anthropogenic determinants of occurrence within an occupancy framework to evaluate habitat suitability of non-protected regions (with sloth bears) in northeastern Karnataka, India. We employed a systematic sampling methodology to yield presence absence data to examine a priori hypotheses of determinants that affected occupancy. These covariates were broadly classified as habitat or anthropogenic factors. Mean number of termite mounds and trees positively influenced sloth bear occupancy, and grazing pressure expounded by mean number of livestock dung affected it negatively. Also, mean percentage of shrub coverage had no impact on bear inhabitance. The best fitting model further predicted habitats in Bukkasagara, Agoli, and Benakal reserved forests to have 38%, 75%, and 88%, respectively, of their sampled grid cells with high occupancies (>0.70) albeit little or no legal protection. We recommend a conservation strategy that includes protection of vegetation stand-structure, maintenance of soil moisture, and enrichment of habitat for the long-term welfare of this species.

Modelling of soldier fly halteres for gyroscopic oscillations

Nature has evolved a beautiful design for small-scale vibratory rategyro in the form of dipteran halteres that detect body rotations via Coriolis acceleration. In most Diptera, including soldier fly, Hermetia illucens, halteres are a pair of special organs, located in the space between the thorax and the abdomen. The halteres along with their connecting joint with the fly's body constitute a mechanism that is used for muscle-actuated oscillations of the halteres along the actuation direction. These oscillations lead to bending vibrations in the sensing direction (out of the haltere's actuation plane) upon any impressed rotation due to the resulting Coriolis force. This induced vibration is sensed by the sensory organs at the base of the haltere in order to determine the rate of rotation. In this study, we evaluate the boundary conditions and the stiffness of the anesthetized halteres along the actuation and the sensing direction. We take several cross-sectional SEM (scanning electron microscope) images of the soldier fly haltere and construct its three dimensional model to get the mass properties. Based on these measurements, we estimate the natural frequency along both actuation and sensing directions, propose a finite element model of the haltere's joint mechanism, and discuss the significance of the haltere's asymmetric cross-section. The estimated natural frequency along the actuation direction is within the range of the haltere's flapping frequency. However, the natural frequency along the sensing direction is roughly double the haltere's flapping frequency that provides a large bandwidth for sensing the rate of rotation to the soldier flies.

An alluaudite Na-2+Fe-2x(2-x)(SO4)(3)(x=0.2) derivative phase as insertion host for lithium battery

A new desodiated derivative compound, Na0.89Fe1.8(SO4)(3), was prepared by the chemical oxidation of alluaudite Na2.4Fe1.8(SO4)(3) Phase using NOBF4 as oxidant. The structure and valency of Fe were characterized by X-ray diffraction (XRD) and Fe-57 Mossbauer spectroscopy. Intercalation behavior of lithium ions in the structure of Na0.89Fe1.8(SO4)(3) was gauged by electrochemical analyses and ex-situ X-ray diffraction. A high capacity of 110 mAh g(-1) at 0.1 C was obtained with a good rate kinetics within a range of 0.1-10 C(1 C = 118 mAh g-1) involving a high Fe3+/Fe2+ redox potential of 3.75 V (vs. Li/Li+). These results confirmed that the Na2.4-delta Fe1.8(SO4)(3) framework was stable even after oxidation and forms a new competitive cathode for the reversible intercalation of lithium ions. (C) 2014 Elsevier B.V. All rights reserved.

The beta hairpin structure within ribosomal protein S5 mediates interplay between domains II and IV and regulates HCV IRES function

Translation initiation in Hepatitis C Virus (HCV) is mediated by Internal Ribosome Entry Site (IRES), which is independent of cap-structure and uses a limited number of canonical initiation factors. During translation initiation IRES-40S complex formation depends on high affinity interaction of IRES with ribosomal proteins. Earlier, it has been shown that ribosomal protein S5 (RPS5) interacts with HCV IRES. Here, we have extensively characterized the HCV IRES-RPS5 interaction and demonstrated its role in IRES function. Computational modelling and RNA-protein interaction studies demonstrated that the beta hairpin structure within RPS5 is critically required for the binding with domains II and IV. Mutations disrupting IRES-RPS5 interaction drastically reduced the 80S complex formation and the corresponding IRES activity. Computational analysis and UV cross-linking experiments using various IRES-mutants revealed interplay between domains II and IV mediated by RPS5. In addition, present study demonstrated that RPS5 interaction is unique to HCV IRES and is not involved in 40S-3 ` UTR interaction. Further, partial silencing of RPS5 resulted in preferential inhibition of HCV RNA translation. However, global translation was marginally affected by partial silencing of RPS5. Taken together, results provide novel molecular insights into IRES-RPS5 interaction and unravel its functional significance in mediating internal initiation of translation.

SInCRe-structural interactome computational resource for Mycobacterium tuberculosis

We have developed an integrated database for Mycobacterium tuberculosis H37Rv (Mtb) that collates information on protein sequences, domain assignments, functional annotation and 3D structural information along with protein-protein and protein-small molecule interactions. SInCRe (Structural Interactome Computational Resource) is developed out of CamBan (Cambridge and Bangalore) collaboration. The motivation for development of this database is to provide an integrated platform to allow easily access and interpretation of data and results obtained by all the groups in CamBan in the field of Mtb informatics. In-house algorithms and databases developed independently by various academic groups in CamBan are used to generate Mtb-specific datasets and are integrated in this database to provide a structural dimension to studies on tuberculosis. The SInCRe database readily provides information on identification of functional domains, genome-scale modelling of structures of Mtb proteins and characterization of the small-molecule binding sites within Mtb. The resource also provides structure-based function annotation, information on small-molecule binders including FDA (Food and Drug Administration)-approved drugs, protein-protein interactions (PPIs) and natural compounds that bind to pathogen proteins potentially and result in weakening or elimination of host-pathogen protein-protein interactions. Together they provide prerequisites for identification of off-target binding.